Electrostatically-Driven Granular Media: Phase Transitions and Coarsening 
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We report the experimental and theoretical study of electrostatically driven granular material. We 
show that the charged granular medium undergoes a hysteretic first order phase transition from 
the immobile condensed state (granular solid) to a fluidized dilated state (granular gas) with a 
changing applied electric field. In addition we observe a spontaneous precipitation of dense clusters 
from the gas phase and subsequent coarsening - coagulation of these clusters. Molecular dynamics 
simulations shows qualitative agreement with experimental results. 

PACS: 47.54.+r, 47.35.+i,46.10.+z,83.70.Fn 



Despite extensive study over the preceding decade, a 
fundamental understanding of the dynamics of granular 
materials still poses a challenge for physicists and engi- 
neers ||^,|[ . Peculiar properties of granular materials can 
be attributed to strong contact interactions and inelastic 
collisions between grains |^,|| . Fascinating collective be- 
havior appears when small particles acquire an electric 
charge and respond to competing long-range electromag- 
netic and short range contact forces. 

The electrostatic excitation of granular media offers 
unique new opportunities compared to traditional vibra- 
tion techniques which have been developed to explore 
granular dynamics It enables one to deal with 

extremely fine powders which are not easily controlled 
by mechanical methods. Fine particles are more sen- 
sitive to electrostatic forces which arise through parti- 
cle friction or space charges in the particle environment. 
Their large surface to volume ratio amplifies the effect 
of water or other surfactants. These effects intervene 
in the dynamics causing agglomeration, charging, etc., 
making mechanical experiments uncontrollable. Electro- 
static driving makes use of these bulk forces, and allows 
not only for removal of these "side effects," but also for 
control by long-range electric forces. 

In this Letter we report experimental and theoretical 
study of electrostatically driven granular material. It is 
shown that the charged granular medium undergoes a 
first order phase transition from the immobile condensed 
state (granular solid) to a fluidized dilated state (granular 
gas) with a changing applied electric field. A spontaneous 
precipitation of dense clusters from the gas phase and 
subsequent coarsening - coagulation of these clusters is 
observed in a certain region of the electric field values. 
We find that the rate of coarsening is controlled by the 
amplitude and frequency of the applied electric field. We 
have also performed molecular dynamics simulations of 
electrostatically driven particles. These simulations show 
qualitative agreement with experiments. 

The experimental cell is shown in Fig. |l|. Conducting 
particles are placed between the plates of a large capaci- 
tor which is energized by a constant or alternating elec- 



tric field. To provide optical access to the cell, the upper 
conducting plate is made transparent. We used 4x6 cm 
capacitor plates with a spacing of 1.5 mm and 35 /im 
copper powder. The field amplitude varied from to 10 
kV/cm and the frequencies on the interval of to 250 Hz. 
The number of particles in the cell was about 10^. The 
experiments were performed both in atmospheric pres- 
sure and in vacuum ( 5 x 10~^ Torr). 

When the conducting particles are in contact with 
the capacitor plate they acquire a surface charge. As 
the magnitude of the electric field in the capacitor ex- 
ceeds the critical value Ei the resulting (upward) elec- 
tric force overcomes the gravitational force mg (to is the 
mass of the particle, g is the acceleration due to grav- 
ity) and pushes the charged particles upward. When 
the grains hit the upper plate, they deposit their charge 
and fall back. Applying an alternating electric field 
E = Eq sin(27r/i), and adjusting its frequency / one can 
control the particle elevation by effectively turning them 
back before they collide with the upper plate. 

The phase diagram is shown in Fig. ^. We have found 
that at amplitudes of the electric field above a second 
threshold value, E2 > Ei, the granular medium forms a 
uniform gas- like phase (granular gas). This second field 
E2 is 50-70% larger then Ei in nearly the whole range of 
the parameters used. In the field interval Ei < E < E2, 
a remarkable phenomenon analogous to coalescence dy- 
namics in systems exhibiting first order phase transi- 
tions was observed. Upon decreasing the field be- 
low E2, the gas phase loses its stability and small clus- 
ters of immobile particles surrounded by the granular 
gas form. These clusters then evolve via coarsening dy- 
namics: small clusters disappear and large clusters grow. 
Eventually the clusters assume almost perfect circular 
form (Figs ||a-c). In the process of coarsening the electric 
current also changes in time due to decrease of particles 
in the gas phase. A close-up image of one of the clusters 
is shown in Fig. |3^. Surprisingly, we can even see a larger 
particle near the center of the cluster which served as the 
nucleation point at the early stages of the aggregation. 
Due to its increased mass the large particle is the first to 
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become immobile in the lower field. After a very large 
time {t « 30000 sec) a single cluster containing about 
10^ grains survives. At the final stage a dynamic equi- 
librium between the granular solid and the surrounding 
gas persists - not all the particles join the last cluster. 

For the cell under atmospheric pressure (open cell) we 
found that both fields Ei and E2 grow as a function of 
frequency for large / and show non-monotonous behav- 
ior for / « 12Hz. This indicates a characteristic time of 
the order 100 msec. We suggest that cohesion may be 
responsible for this relatively large time. Indeed, due to 
the humidity of the air a surface coating should exist, 
requiring a characteristic time r for the grain to detach 
from the capacitor plate. In order to reduce the cohesion 
we evacuated the cell to 5 x 10^^ Torr. As demonstrated 
from Fig. |^, the frequency dependence is indeed substan- 
tially reduced and becomes almost flat. Small oscillations 
in the dependence for low frequency are probably due to 
a residual coating on the particles which does not com- 
pletely evaporate in vacuum. 

We have measured the number of clusters N and the 
averaged radius of the clusters (i?) as a function of time t 
(see Fig. ^ using images taken in time-lapse by a digital 
camera and then processed on a computer. The values 
of (R) for both open and evacuated cells are consistent 
with a power law {R) ^ \fi. However, there is a sub- 
stantial difference in the behavior of TV. In open cells 
(no vacuum) we observed slow coarsening N ~ 
whereas the coarsening for the evacuated cell is much 
faster, iV ^ Xjt. Remarkably, for the evacuated cell the 
exponents for N and (i?) coincide with the exponents 
for two-dimensional bi-stable systems in the interface- 
controlled regime of Ostwald Ripening |^,|l^ . We specu- 
late that the slow-down in coarsening for open cells can 
be a result of cohesion between the particles and with 
the capacitor plate. This cohesion reduces the mobility 
of particles at the edges of the clusters, resulting in "pin- 
ning" of the front between the granular solid and gas. 

Let us compare the forces exerted on an isolated par- 
ticle and on a particle which is held within a cluster at 
the same field. The force between an isolated sphere and 
the capacitor plate in contact can be found as a limit of 
the problem of two spheres when the radius of one sphere 
goes to infinity. This problem had been contemplated by 
several outstanding physicists including Kelvin, Poisson, 
Liouville and Kirchhoff Building on the technique 

of Ref. ||ll|], we arrive at the force in question: 

Fe = ca'E^ (1) 

where a is the radius of the sphere and E is the field in 
the capacitor. The constant c = ((3) -I- 1/6 ~ 1.36 comes 
from summing of infinite series derived in . The elec- 
tric force Fg has to counterbalance the gravitational force 
G = A/Snpga^, where p is the density of the material. 
Comparing and G we find the first critical field Ei: 



Our theory indicates no frequency dependence of Ei . For 
the parameters of our experiments we obtain Ei = 2.05 
kV/cm. The experimental critical field is Ei « 2.4 
kV/cm. This discrepancy seems reasonable since we did 
not take into account additional molecular and/or con- 
tact forces which will increase the critical field. 

If the spheres are close to each other, the surface 
charge will redistribute: each sphere in the cluster ac- 
quires smaller charge than that of an individual sphere 
due to a screening of the field by its neighbors. The exact 
derivation of the force acting on the particle within the 
cluster is not available at present. The upper bound can 
be obtained by replacing the square lattice of spheres 
with radius a by a lattice of squares with area 4a^. 
Since the charge density of the corresponding flat layer 
is fj = F/47r, the total electric force on the square is 
F2 Aa^(jE/2. The ratio of the fields to hft indi- 
vidual particle Ei and the particle in the square lat- 
tice E2 is E2IE1 = \/27r • 1.36 ^ 2.92. For closed- 
packed hexagonal lattice one obtains a slightly higher 
value E2/E1 3.14. As one sees, the ratio of the critical 
fields is independent of particle size and the density and 
the frequency /. This is consistent with the experimental 
findings, see Fig. ||, although it exceeds the experimental 
value by a factor of 2. A more accurate account for the 
surface shape will improve the ratio. 

We have also performed molecular dynamics simula- 
tions of conducting spheres in applied electric field. Sev- 
eral simplifications have been implemented: the sphere 
polarization was neglected; the charge was assumed to 
be in the center of the sphere; all collisions were assumed 
to be inelastic. These simplifications are justified by the 
fact that the particle-particle collisions are rather rare 
in contrast to deep vibrated granular layers studied in 
Ref. Interactions between the spheres are treated 

as being between point charges if the distance is smaller 
than the plate spacing. For larger distances we assumed 
no interaction due to screening of the far charge field by 
other particles (compare with Debye screening |8|). 

In our molecular dynamics simulations we imple- 
mented explicitly that the charge acquired by each sphere 
from the conducting plate decreases gradually with the 
number of nearest neighbors: it is maximum for an iso- 
lated sphere and decreases by factor of two for the sphere 
inside the cluster. We simulated 35 /im diameter 3000 
copper spheres in the domain 7x7 mm (200 x 200 particle 
diameters) with periodic boundary conditions. The ca- 
pacitor spacing was 1.5 mm, as it was in the experiment. 
Although the quantitative correspondence between the 
molecular dynamics simulations and the experiment was 
not planned (the number of particles and the size of the 
system are too small ) , the simulations turns out to be in 
a qualitative agreement with the experiment. We have 
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obtained a gas phase for high field levels and spontaneous 
formation of clusters and coarsening for lower filed level, 
see Fig. ||. However, detailed simulations with much 
large number of particles and more realistic account for 
sphere polarization effects are necessary to achieve quan- 
titative agreement with experiment. 

We have reported a hysteretic first-order phase tran- 
sition accompanied by coarsening behavior in the 
electrostatically-driven conducting granular medium. 
The origin of coarsening dynamics and hysteresis is due 
to a screening of the electric field in dense particle clus- 
ters. Our results show surprisingly high sensitivity of 
the phase boundaries to surface coating of the grains due 
to humidity. Simplified molecular dynamics simulations 
with conducting particles demonstrate qualitative agree- 
ment with the experiment: existence of two critical fields 
and coarsening. 
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FIG. 1. The electrostatic cell 

FIG. 2. The phase diagram in /, _Eo plane. Lines with 
diamonds show i5i,2 in open cell, dashed (dash-dotted) lines 
show E2 (El) for the evacuated cell. In gas only part clusters 
do not form, in coarsening part clusters and gas coexist. 

FIG. 3. (a) Close-up image of a single cluster. The larger 
particle near the center (arrow) of the cluster served as the 
nucleation point. Top view of the cell for (b) t — 10 sec, and 
(c) t = 10910 sec for / = 40 Hz and E = 5.25 kV/cm. 

FIG. 4. (R) vs t^/^ for open cell (o), evacuated cell (•) 
and 1/A'' for open cell (o). The conditions were E = 5.25 
kV/cm and / = 40 Hz for open cell and E = A kV/cm and 
f — 8 Hz for evacuated cell. Dashed lines show linear fit to 
these curves. Inset: N vs t for open cell (o). Dashed lines 
show power-law fits ~ and N ~ l/t correspondingly. 

FIG. 5. Top view of simulation domain. Size of the do- 
main is 200 X 200 particle diameters, number of particles is 
3000, value of the applied field is -E = l.l_Bi. The configu- 
rations are shown for t — Isec (a) and t = lOsec (b). Black 
bullets correspond to immobile particles in the clusters, open 
circles correspond to fiying particles. 
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